Study of Photonuclear Interaction of Muons in rock 
with the MACRO experiment. 
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We present first results about the measurement of the charged hadrons production by atmospheric muons in the 
rock above MACRO. A comparison between the measured rate with the Monte Carlo expectation is presented. 



1. Introduction 

The inelastic muon-nucleus interaction was 
studied at accelerators, in the range of large Q^ 
transfer (Q^ > IGeV^) mainly to measure the 
nucleon structure functions (deep inelastic scat- 
tering experiments). However, the bulk of in- 
teractions are characterized by low Q^ (Q^ < 
O.lGeV^) and can be described with the exchange 
of a quasi-real photon between the muon and 
the nucleon and they are often referred to as 
photo- nuclear interaction of muons. Recently, 
it has been stressed that nuclear interactions of 
muons are an important source of background for 
many underground experiments. A comparison 
between experimental measurements about the 
hadrons production by muons and Monte Carlo 
(MC) simulations is mandatory to test the relia- 
bility of theoretical models. In fact, uncertainties 
do exist on the cross section calculation for the 
process, mainly due to the extrapolation of photo- 
nuclear cross section at high energies and in the 
simulation of the hadronic final state, because the 
smallness of the Q^ does not allow the use of per- 
turbative QCD. The aim of the present work is to 
compare real data collected by the large area un- 
derground experiment MACRO ||l|, taking data 
at Gran Sasso laboratory, with GEANT 3.21 and 
with FLUKA97, the latter based on the cross sec- 
tion calculated by Bezrukov and Bugaev||^ and on 
DPM modelW for the sampling of the final state. 



*For the complete list of the Collaboration see the paper 
" Relevance of the hadronic interaction model in the inter- 
pretation of multiple muon data as detected with MACRO 
experiment" by O. Palamara at these proceedings 



2. The Detector 

A description of MACRO is given elsewhere Q|. 
We remind here that MACRO operates at an av- 
erage depth of 3800 hg/cm^, where the average 
residual energy of muons is about 300 GeV. The 
horizontal area of MACRO is about 1000 m^ and 
charged tracks can be reconstructed by means of 
a streamer tube system, in two projective views, 
with a space point accuracy of ~1 cm. The up- 
per part of MACRO is in practice a thin detector, 
imposing a low threshold on secondary particles. 
The lower part, thanks to the rock absorbers, al- 
lows to stop particles and to measure their range 
up to few hundreds of MeV. We consider events 
in which the muon interacts in the rock above 
the apparatus and both the muon and at least 
one charged hadron are observed in the horizon- 
tal planes of the tracking detector. An example 
is shown in Fig.l, in which a muon enters the 
apparatus from above, along with at least one 
charged hadron. Muon tagging is performed us- 
ing the standard MACRO muon tracking while 
the additional tracks have been recognized with 
a specialized algorithm, looking for shorter tracks 
pointing towards the main track around a com- 
mon vertex region contained in the rock above the 
detector. In practice we select charged hadrons 
with a minimum kinetic energy around 150 MeV. 

3. The physical process and its simulation 

We have chosen FLUKA as a main reference 
to generate the underground muon events, taking 
into account the following steps. The direction of 
incident muons is sampled according to the local 
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For a more detailed discussion of the differences 
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Figure 1. A candidate event detected in MACRO. 

angular distribution measured by MACRO and 
properly unfolded. As far as muon residual en- 
ergy at a depth h (in km w.e.) is concerned, we 
chose to sample it according tog 

dN(E^,h) 
d(Ej 



= Ke 



-bh(7-l) 
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where 7=3.5,b=0.4(km w.e.)^^ and e w540 GeV. 
Given the direction, the slant depth of rock h 
crossed by the muon can be obtained from the 
map of the mountain overburden. In the simu- 
lation, muons are allowed to interact in a 13 m 
thick layer of rock positioned all around the ex- 
perimental hall. This thickness corresponds to 
about 35 Xint for hadrons, and this is enough 
to fully contain hadronic showers. If a photo- 
nuclear interaction occurs in the region of rock 
described above, the muon and secondary parti- 
cles are transported through the rock, along with 
e.m. and hadronic showers possibly produced. If 
both the muon and at least one additional parti- 
cle reach the tunnel, the event is stored and a full 
simulation of the apparatus response is performed 
using the GEANT 3.21 packageg]. The simu- 
lated data are treated with the the same analy- 
sis instruments of real data. In order to make 
a comparison with a different model, we have 
repeated the event generation in the rock with 
GEANT. One of the most important differences 
is in the total cross section, as shown in Fig. |[ 
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Figure 2. The total muon photonuclear cross sec- 
tion in the considered models. 

between FLUKA and GEANT-GHEISHA for this 
process, see 0. 

4. Data Analysis 

The main difficulty in our analysis is to achieve 
the necessary rejection factor against the physical 
background. Such a background is largely domi- 
nated by two processes: 1) the e.m. interactions 
of muons in the rock (bremsstrahlung and pair 
production, which have a much larger cross sec- 
tion than the photonuclear reaction); 2) multi- 
ple muon events in which one of the muons stops 
inside the detector. An example of e.m. back- 
ground event is shown in Fig. H. This back- 
ground has been extensively studied with our sim- 
ulation tools. In the case of the e.m. interactions, 
the mean angular separation between the muon 
and the additional tracks is less than the corre- 
sponding separation observed in photonuclear in- 
teractions. Besides that, e.m. events often show 
very large clusters of fired tubes near the muon 
track. These features are used to achieve a re- 
jection factor against e.m. events at the level of 
4.5- 10~^, maintaining the recognition efficiency 
for hadronic events at the level of 55%. In or- 
der to reject the muon bundle backgrounds, ad- 
ditional cuts based on the parallelism of tracks 
have to be considered. We have studied the events 
generated with the HEM AS code||] using for the 
primary cosmic ray spectrum and mass composi- 
tion the results from the best fit reproducing the 
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Figure 3. Example of background event: e.m. in- 
teraction of muon in the rock. 
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Figure 4. Charged hadrons multiplicity distribu- 
tion 

MACRO data themselves Q The rejection factor 
achieved for muon bundle background at the level 
of 1.5-10"'*, at the price of a slight reduction of 
the selection efficiency, which is now at 47%. 

5. Results 

We have analysed a data sample corresponding 
to about 11000 hours of full running of the de- 
tector. With the above selection criteria, we have 
found 1938 candidate events over a total sample 
of 9544318 muon events. From our knowledge of 
the background, we expect a contamination by 11 
events from the e.m. interactions in the rock, and 



by 107 events muon bundles surviving the cuts. 
We can express the results in terms of the ra- 
tio Rfj,+h of the selected /x-|-hadrons events to the 
number of muon events in the same live time. We 
then compare the experimental results to the MC 
prediction having used the same selection criteria. 
After the subtraction of the background, we find 
for Rfj^+h in real data and in the MC simulations 
the following results: 



jdDATA 
nFLUKA 



(1.91±0.05,ta±0.03,j,,)-10- 

(1.89±0.16,ta±0.02,j,,)-10- 



- R^+h""^ = (1-31 ± 0.14,ta ± 0.02,,,) . 10-5. 
The systematic error on the experimental data 
is due to the uncertainties on background sub- 
traction, while in the simulation is dominated by 
the uncertainties on the muon energy spectrum. 
Fig. 4 shows the charged hadron multiplicty dis- 
tribution. A good agreement between the real 
data and FLUKA is obtained, while the GEANT 
distribution is flatter, reflecting the Q^ spectrum 
harder than FLUKA W. From this preliminary 
measurement, we can conclude that the MACRO 
experiment can measure of the charged hadron 
(Efei„>150 MeV) production in the rock at the 
desired level of accuracy. The FLUKA predic- 
tions, based on the Bezrukov and Bugaev model, 
are in very good agreement with data, while the 
GEANT-GHEISHA model gives absolute predic- 
tions lower by an order of magnitude. 
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